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ABSTRACT 

In  9  medium  in  Which  noise  and  signal  propagate  with  different  velo¬ 
cities,  a  detection  scheme  making  use  of  the  cross-correlation  between  the 
received  waveforms  at  two  different  locations  offers  a  number  of  distinct 
advantages  over  one  making  use  of  ordinary  autocorrelation  of  the  received 
waveform.  These  advantages  arise  from  the  resultant  relative  shift  between 
the  noise  and  signal  autocorrelation  functions  allowing  for  separation  of 
noise  and  signal  when  the  noise  characteristics  or  statistics  are  unknown, 
and  even  where  the  processes  are  non -stationary.  Ramifications  of  the  cor¬ 
relation  shift  are  an  Improvement  in  signal-to-noise  ratio  (even  when  noise 
characteristics  are  known),  and  capability  for  measurement  of  velocity,  sig¬ 
nal  source  direction,  and  signal  and  noise  and  channel  characteristics. 
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PREFACE 

With  the  current  interest  in  detection  of  underground  nuclear  explosions, 
this  report  has  been  prepared  to  point  out  that  effective  use  can  be  made  of 
the  dual -velocity  properties  of  the  medium.  The  report  is  a  modification  of 
a  memorandum  originally  prepared  in  1960  while  at  the  Applied  Science  Divi¬ 
sion  of  Melpar,  Xne.,  Watertown,  Massachusetts.  X  would  like  to  thank  James 
Lee,  presently  of  Raytheon  Company,  for  pointing  out  that  seismic  signals 
propagate  with  different  velocities  and  suggesting  that  two  receivers  might 
be  used  to  exploit  this  property. 
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I.  INTRODUCTION 

Usual  detection  of  a  signal  in  noise  involves  cross-correlation  when  the 
expected  received  signal  waveform  is  known,  and  auto-correlation  when  it  is 
unknown,  all  at  a  single  receiver.  There  are  environments  in  which  signal 
and  noise  propagate  with  different  velocities.  A  dual  receiver  detection 
scheme  in  Which  the  received  waveforms  at  two  different  locations  are  cross- 
correlated  offers  a  number  of  distinct  advantages  over  ordinary  auto-correla- 
tion  in  detection  of  signals  having  traversed  a  multi -velocity  medium.  This 
Note  describes  some  of  these  advantages,  particularly  in  terms  of  a  duo- 
velocity  medium  in  Which  noise  and  signal  propagate  at  two  different  velo¬ 
cities,  and  where  the  problem  is  one  of  passive  detection,  communications,  or 

1  *8 

measurement  of  signal  or  channel  characteristics.  Previous  examination  of 
dual  receivers  has  not  focussed  on  a  velocity-differential  environment,  and 
phased-array  studies  have  used  diversity  combining  drawing  upon  addition  of 
many  receiver  outputs  rather  than  correlation  of  two. 

Aa  is  well  known  the  seismic  channel  supports  multi -velocity  propagation 

*  9-11 

modes  ’  ,  and  this  could  be  of  special  interest  in  detection  of  under¬ 

ground  nuclear  explosions  or  earthquakes  or  in  differentiating  the  two  (even 
for  simultaneous  occurrence)  ,  in  underground  communications,  and  in 

signal  probing  to  study  the  Earth  structure.  Interest  may  center  on  differ¬ 
entiating  between  two  aperiodic  or  two  stochastic  waveforms,  or  on  detecting 
a  signal  imbedded  in  the  natural  mieroseismic  noise  of  the  Earth. 

Perhaps  underwater  sound  propagation  and  electromagnetic  waves  in  the 
ionospheric  r  .annel  exhibit  similar  anomalies  in  velocity  of  propagation. 

This  paper  treats  the  properties  of  the  dual- receiver  detection  scheme,  with¬ 
out  further  reference  to  the  apacific  propagation  media  to  which  it  can  be 
applied. 

V . . . 

As  described  in  some  of  the  references,  the  principal  types  of  propagation 
modes  (having  different  velocities)  in  the  seismic  channel  are  the  compres¬ 
sions!  (P)  waves  and  the  shear  (B)  waves  in  the  body,  and  the  bov®  anTi 
Kay t nigh  surface  waves,  the'  ftaylelgh  accounting  for  much  of  the  micronetnmic 
noise,  Boms  of  the  discussion  of  this  paper,  e.g.,  that  dealing  with  direc¬ 
tion  finding,  assumes  e  somewhat  idealised  channel  mode]  in  which  the  propa¬ 
gation  characteristics  fcr  each  mode  are  assumed  essentially  invariant  through¬ 
out  the  propagation  path.  However,  the  main  discussion  concerning  detection 
assumes  only  that  the  propagation  characteristics  remain  invariant  in  the 
region  between  the  two  component  receivers  of  the  dual  receiver,  and  the 
velocity  components  colinear  with  the  axis  of  the  dual -receiver  dipole  ore  the 
ones  of  interest.  Actually,  even  if  this  latter  invariance  did  not  exist,  the 
mnlor  benefits  of  the  relative  correlation  shift,  would  still  arise. 
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The  problem  of  interest  is  one  in  which  we  wish  to  detect,  at  the  re¬ 
ceiver,  the  presence  of  signal  (and  noise),  as  compared  with  just  noise,  vhcr 
the  noise  is  microselsmic  or  any  undesired  signal.  The  resulting  overall 
cross-correlation  function  at  the  output  of  the  dual  receiver  consists  of  the 
signal  auto-correlation  function,  the  noise  auto-correlation  function  dis¬ 
placed  from  the  former,  and  finite-time  correlation  noise.  The  advantages 
derived  with  the  dual-receiver  all  arise  from  the  relative  shifting  of  the 
signal  and  noise  correlation  functions.  The  discussion  will  treat  examina¬ 
tion  of  the  received  waveform  at  a  specified  time,  at  which  time  the  overall 
correlation  function  is  sampled  at  the  peak  of  the  signal  auto-correlation 
function.  In  contrast  to  ordinary  auto- correlation ,  the  noise  car  be  separa¬ 
ted  from  the  signal  term  when  the  noise  statistics  are  not  adequately  known. 
Even  when  the  noise  statistics  are  known,  there  can  be  a  4  db  improvement  in 
effective  signal-to-noise  ratio  over  that  obtained  from  auto-correlation,  in 
a  passive  system  where  it  is  not  possible  to  simply  increase  the  transmitter 
powor,  this  saving  may  be  the  final  factor  which  influences  whether  or  not-, 
the  required  detection  reliability  is  achieved. 

II.  SHIFTED  CORRELATION  FUNCTIONS 

We  treat  an  Idealized  model  to  emphasize  the  method.  The  dual  receiver 
has  the  configuration  shown  in  Figure  1.  The  waveform  is  received  at  two 
points,  and  the  channel  between  the  two  receivers  is  treated  as  a  bi -modal 
delay  line  ,  having  a  delay  of  Tg  for  the  signal  and  T  for  the  noise. 

■*  . . 

This  paper  hr  been  concerned  with  general  conclusions,  rather  then  with 
fine  detail.  If  there  is  channel  distortion  between  the  two  receivers, 
then  of  course  the  expressions  for  uz  and  of  must  be  modified.  The  inter¬ 
vening  channel  could  then  be  represented  byztwo  delay  lines,  one  each  for 
signal  end  noise,  and  each  followed  by  filters  having  impulse  responses 
respectively  of  ha(t)  and  hn(t),  (These  two  impulse  responses  could  be  tin 
same.) 

We  then  have 

ux  "  J  V0)VT  "  °*d0  f  J  Rns*°*h.M*T  ‘  c,*dd 

v  f  R.-tcOK.t*  -  L  -  o)da  I  R  (c)h  (t  -  L  -  c)do. 
j  »n  n  j  n  n 

Tin.  expressions  derived  in  the  Appendix  for  variance  and  hence  (fl/N)  are 
?  \m  I  arly  altered.  " 
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THE  (INTERVENING)  CHANNEL 


Figure  1,  Dual  Receiver 
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(The*®  delay*  reflect  the  different  velocities  for  signal  end  noise) .  The 
waveform  at  receiver  1  i*  x(t)  ~  s(t)  +  n(t)  ,  where  the  component  terms  de¬ 
note  signal  and  noise  respectively.  At  receiver  2  the  waveform  is 
y(t)  =  s(t  -  Tfl)  +  n(t  -  TR) .  The  detector  consists  of  an  artificial  delay, 
and  a  finite-time  correlator.  For  simplicity  in  this  discussion  we  examine 
low-pass  correlation  (for  video  waveforms)  ,  and  the  correlator  then  consists 
of  a  multiplier  followed  by  a  finite-time  integrator.  The  output  of  the 
correlator  is  sampled  for  decision  at  the  time  at  which  the  signal,  if 
present,  has  effectively  been  completely  received  at  the  second  receiver. 

The  correlator  output  can  be  a  random  variable16’17,  especially  in  the 

presence  of  noise,  and  this  output  consists  of  the  average  component  (which 

is  the  one  of  interest)  and  a  fluctuation  component.  (See  Appendix) .  The 

average  output  is  then  u„  =  R.(t)  +  R  „(x)  +  R  (x  +  T-  T)+R(x  +  T  -T), 

3  zb  ns  sn  s  n  n  an’ 

where  Rfi  and  Rn  are  the  signal  and  noise  auto-correlation  functions,  respec¬ 
tively,  and  the  other  two  terms  are  the  signal  and  noise  cross-correlation 
functions.  When  *(t)  and  n(t)  are  uncorrelated  we  obtain  u.  -  R_(x)  t  r  (t  -  L) 

where  L  ■  T  -  T  »  and  in  particular  at  x  ~  0  this  becomes  u_  ”  R  (0)  t  R  (L)  , 
n  m  z  a  n 

We  examine  the  function  for  x  =  0,  and  when  L  is  large  enough,  and  hence 

Rr(L)  is  small  enough,  we  effectively  sample  only  the  signal  auto-correlation 

function.  In  Figure  2  we  show  u  consisting  of  the  signal  and  the  displaced 

z 

noise  auto-correlation  functions,  and  in  Figure  3  we  show  the  same  for  ordin¬ 
ary  auto-correlation  of  the  received  waveform,  corresponding  to  a  zero  dis¬ 
placement,  L  =  0. 

Normally  when  autocorrelation  is  used  the  noise  statistics  are  also  at 
least  partially  xnown,  and  in  particular  Rn(0)  is  known  and  can  be  subtracted 
out  from  the  u_  leaving  just  R.(0) .  Thus  examination  of  u_  -  R_(0)  indicates 
presence  or  absence  of  signal,  depending  upon  whether  or  not  this  is  non-zero 
(considering,  for  the  moment,  the  average  value,  and  ignoring  the  fluctuation 
component).  Thus,  if  we  know  Rn(0)  then  it  can  be  subtracted  out  anyway,  and 
there  is  then  no  direct  advantage  in  obtaining  a  displacement  between  the  two 


Even  when  signal  and  noise  are  correlated  we  get  proper  indication  of 
whether  signal  is  present  or  not.  The  function  uz  then  consists  of  all 
four  terms  indicated  above,  but  as  also  show  above  only  the  first  two 
terms  are  significant  at  x  =  0  when  L  is  larqo  enough.  Thus  at  x  “  0 
the  significant  terms  are  R„(0)  +  RnB(0)  »  and  both  are  present  when  signal 
is  present  and  both  are  absent  when  signal  is  absent. 
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correlation  functions.  (There  still  is  however,  an  indirect  advantage  re¬ 
sulting  in  Improved  signal -to-noise  ratio,  and  this  i»  cited  in  the  next  nub- 
section.)  However,  when  the  noise  statist icr  are  not  known  and  when  Rn(0) 
cannot  be  subtracted  out,  there  great  merit  in  being  able  to  achieve  a 

relative  displacement  in  the  correlation  functions. 

The  correlation  shift  is  effective  with  either  stochastic  or  determinis¬ 
tic  waveforms.  It  is  of  particular  value  for  non-stationary  processes,  for 
which  the  noise  term  could  not  ordinarily  be  subtracted  out,  the  statistics 
themselves  fluctuating.  Although  we  have  described  the  problem  for  known 
arrival  times,  the  correlation  shift  is  valuable  also  for  continuous  scanning 
of  the  waveform,  since  a  new  correlation  function  appears  when  desired  signal 
is  present. 

III.  SIGNAL-TO -NOISE  IMPROVEMENT 

The  probability  of  correct  decision  at  the  output,  of  the  correlator  de¬ 
pends,  of  course,  upon  the  effective  signal -to-noise  ratio  at  that  point.  If 
we  know  the  noise  statistics  and  can  therefore  subtract  out  Rn ( L) ,  the  signal- 
to-noise  ratio  at  the  correlator  output  is 


2 

where  0Z  is  the  variance  of  the  correlator  output.  The  (S/N)q  is  developed 
in  the  Appr  Jix  for  increasing  amounts  of  assumptions  and  specializations. 
For  stochastic  signal  and  noise  which  are  Gaussian  and  uncorrelated,  and 
where  the  noise  is  white  and  band-limited  at  W,  and  (S/N),  is  small,  we  ob¬ 
tain 

<S/N)q  -  2  TW  <S/N)2  [l  ♦  a-fV,rt^T~]  1  * 

where  the  subscripts  i  and  o  refer  respectively  to  the  correlator  input  and 
output.  Maximum  (S/N)0  occurs  for  L  •  ( 3/8W) ( 257 . 5/270) ,  and  it  is  4  db 
better  than  the  value  for  L  ■  0,  the  latter  corresponding  to  ordinary  auto¬ 
correlation.  The  maximum  is  essentially  satisfied  at  i.  ( 3/0W) ,  and  lor 
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known  velocity  differential  thio  can  be  achieved  simply  by  setting  the 
separation  distance  of  the  receivers  as 


where  v_  and  v_  are  velocities  of  propagation. 

m  n 

We  note  that  for  L  equal  to  an  integer  multiple  of  (1/4W)  or  for 

L  »  (1  /VI)  ,  the  improvement  in  signal -to-noise  ratio  over  the  value  for 

L  -  0  is  3  db.  We  might  just  remark  that  for  noise  band-limited  at  w,  R  (L) 

n 

has  zeros  for  L  an  integer  multiple  of  (1/2W) .  Then  Rn(L)  disappears  of 
itself  (at  x  =  0) ,  even  if  we  don't  know  the  actual  values  for  the  noise 
correlation  function.  However,  we  then  achieve  only  the  3  db  improvement  in 
signal -to-noise  ratio. 

IV.  CORRELATION,  VELOCITY,  AND  DIRECTIONAL  MEASUREMENTS 

Although  our  major  interest  has  focussed  on  signal  detection ,  various 
related  measurements  can  be  made  with  the  receiver.  For  example,  even  if 
noise  and  signal  are  correlated,  then  knowledge  of  one  correlation  function 
can  permit  determination  of  the  others.  If  Rn(x)  is  known,  then  from  samples 

of  the  function  U_  =  R  (x)  +  R.(x)  *  R_(x  -  L)  +  R(x  -  L)  around  x  ~  L,  we 

z  s  no  an  n 

can  determine  R  (x)  ,  since  the  first  two  terms  are  negligible  at  t  L  when 

Bn 

L  is  large  enough.  From  this,  Rns(x)  is  known,  and  then  from  examination  of 
u^  in  the  vicinity  of  x  =  0,  we  can  determine  Rs(x).  For  stationary  stochas¬ 
tic  noise,  Rn(x)  could  readily  be  determined  from  measurements  prior  to  sig¬ 
nal  reception 

Signal  probing  can  be  used  for  determination  of  the  velocity-differential 
characteristics  of  the  medium.  For  example,  an  active  system  can  generate  a 
strong  signal  which  may  propagate  over  different  modes,  and  the  velocity 
differential  determined,  W*  might  note  that  two  sources  generating  slqnals 
Which  propagate  at  the  same  speed  but  which  impinge  upon  the  receiver  dipole 
from  different  directions  give  rise  to  an  apparent  velocity  differential. 

For  the  moment  we  digress  and  comment  on  an  application  where  a  dual 
receiver  has  of  course  been  used  for  direction  finding  and  source  location, 
wherein  a  single  source  is  located  as  being  on  a  particular  hyperbola.  In 
our  cases  of  interest  the  distance  between  receivers  is  small,  and  a  signal 
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source  at  a  much  greater  distance  is  located  precisely  in  direction,  since  it 
will  be  lying  essentially  on  the  asymptote  of  one  of  the  aforementioned  hyper¬ 
bolas,  all  asymptotes  bisecting  the  axis  of  the  dual  receiver.  (An  assump¬ 
tion  here  is  straight -line  propagation  and  a  broad  wave  front) .  The  direction 
of  the  signal  source  is  described  by  0,  the  angle  between  the  perpendicular 
to  the  dual  receiver  axis  and  the  line  from  its  center  to  the  source?  by  v, 
the  speed  of  propagation?  D,  the  distance  between  receivers;  and  A,  the  delay 
introduced  in  the  correlation.  Then  D/A  -  v  sin  0,  and  knowing  either  v  or 
0,  the  other  can  then  be  determined.  With  a  triple-receiver,  both  v  and  0 
can  be  determined. 


V.  REMARKS 

Under  certain  circumstances  the  different  velocity  components  could  be 
filtered  out  with  differently  oriented  or  with  different  types  of  sensors. 
However,  When  the  filtering  must  in  essence  be  done  with  the  data  processing 
rather  than  with  the  sensors,  the  detection  procedure  described  has  the 
values  discussed  and  which  are  summarized  below: 

1)  There  is  the  relative  shift  in  signal  and  noise  correlation  functions, 
arising  for  deterministic  and  for  stochastic  signals,  and  particular¬ 
ly  effective  for  signal  detection  when  the  noise  statistics  are  un¬ 
known.  Of  special  interest  is  the  case  where  noise  (and  signal) 
statistics  are  non -stationary ,  and  Rr(0)  could  not  be  subtracted 

out.  Even  aperiodic  noise  and  signal  occuring  at  the  same  time 
(e.g.,  simultaneously  occuring  earthquake  and  nuclear  explosion)  can 
be  separated. 

2)  There  is  the  improvement  in  signal-to-noise  ratio. 

3)  Auto-correlation  and  cross-correlations  functions  can  be  determined 
just  knowing  one. 

4)  There  is  application  to  velocity  determination,  structural  study, 
and  direction  finding. 

The  duo-velocity  discussion  extends  easily  to  the  multi -velocity  case. 
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APPENDIX 

We  derive  the  signal -to-noise  relation  at  the  output  of  the  cross-corre- 
iator  for  the  dual -receiver  detection  scheme.  The  following  discussion  is 
developed  with  increasing  specialization,  and  for  any  actual  problem  the 
degree  of  generality  desired  will  Indicate  how  far  one  should  go  in  using 
these  relations.  Intermediary  steps  are  included  to  allow  one  to  use  the 
expressions  relevant  to  one's  problem.  The  derivation  here  given  is  for  sig¬ 
nal  and  noise  processes  Which  are  random  and  stationary.  An  ensemble  aver¬ 
aging  is  Indicated  by  <  >.  We  will  obtain  the  output  signal -to-noise  for 
t  =  0,  corresponding  to  the  peak  of  the  signal  correlation  function. 


t 

z  ( t  ;t)  =  jp  J  x(t-r-Ta)  y(t)dt. 

t-T 


ub(t)  =  <z ( t f t )  >  =  if  J  <x(t-t-Ts)  y (t)  >  dt 


=  R  (t)  +  R__(t)  +  R  (t -L)  +  R  (t-L)  . 

■  ns  sn  n 

» 

When  n(t)  and  s(t)  are  uncorrelated,  u_(x)  =  R  (x)  +  R  (x-L)r  and  in  particu- 

z  s  n 

lar,  at  x  »  0,  this  is  R„(0)  +  R(L),  which  we  denote  u„. 

T  ’  m  n  z 


[-.<*>]*  “  <z2(t,x)>  -  [uB^)]2 
T  T 

J  J  <x(t-r*T(|)  y(t)x(p-t-Tg)  y((-)>dtdp  -  f ux (x ) I 


T* 


o  o 


ft. 
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At  x  =  0  this  becomes t 
T  T 

k<°,]2  =d2=  \  ]  J  <  [s(t)+n(t)]  [s(t)+n(t-L)]  [s(p)+n(p)]  [s(p)+n(p-L)]  >dtdp-u^ 
T  o  o 

Now  assuming  s(t)  and  n(t)  to  be  Gaussian  with  zero  mean,  thereby  permitting 
use  of  the  relation  <x^x2x3x^>  =  <x^x2><x3x4>  +  <XjX3Xx2x4>  +  <x,x4Xx2x3>, 
and  further  assuiTiing  s(t)  and  n(t)  are  independent  (uncorrelated)  ,  thereby 
allowing  us  to  discard  a  greater  number  of  terms,  and  making  use  of  the 
relation 


T  T  T 

T  J  J  R(Y"P)dYdP  ■  f  J  U  ■  f )  R(a)da,  we  obtain* 

T  o  o  o 

«l  =  i  J  <1  -  T>  (i  -fl  [*.<•>  -  a« 

o  o 

T  T 

+  J  J  (1  -  ip)Rn(a+L)Rn(a-L)da  +  £  J  (1  -  |)Rs(a) [Rn(a-L)  +  Rn(a+L)]da. 
o  o 

Let  us  consider  the  case  in  which  R_(a)  and  R  (a)  -*  0  for  a  >  A  <<  T.  This 

s  n 

permits  dropp^  g  the  (1  -  a/T)  for  the  first  two  integrals.  If  L  >  A,  then 
the  third  integral  vanishes.  We  are  interested  in  retaining  the  last  two 
integrals;  therefore,  let  us  treat  the  case  where  L  is  not  too  large,  say 
L  <  A,  Then  (1  -  a/T)  may  be  dropped  completely.  Keeping  these  modifica¬ 
tion*  in  mind,  we  can  let  the  upper  limit  be  infinity  with  negligible  change 
in  the  result.  Further,  since  the  integrands  are  even,  we  can  express  the 
integrals  over  the  entire  line.  With  these  changes  the  noise  is 

. . . .  2 

When  L  ■=■  0,  corresponding  to  auto-correlation,  ox  becomes 

T 

T  I  *  T5  f  %(a)  *'  Vl,)]3  rta* 
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oo  oo 

al  =  T  J  r2u)<s«  +  t  J  [va)  +  Va)]2  <5° 

-00  -oo L 

oo  oo 

+  i  J  Rn(a+L)  Rn(a-L)da  +  ^  J  Rg(a)  |Rn(a-L)  +  Rn(a+L)l  da. 

-oo  -oo  L  J 

The  Fourier  Transforms  of  Rg(p)  and  Rn(p)  are,  respectively,  G(f) ,  the  signal 
power  spectral  density,  and  N(f) ,  the  noise  power  spectral  density. 


oo  CD 

a2  ■  §  J  02(f)df  +  J  J  N2(f)  [1  +  cos  UrrLf)]  df 
-00  -00 


00 

+  J  J  0(f)  N(f)  [l  +  cos  (2TTLf)]  df. 
-00 


For  signal  and  for  white  noisa  (of  density  NQ)  both  band-limited  at  W j 


+W 

T  J  o2(«)d«  + 

-w 


Khan  0(f)  <<  N0,  tha  first  tarm  is  negligible,* 

In  terms  of  tha  input  signal -to-noiae  ratio,  (S/N)i  Ra(0)/2WNQ,  the 
output  signal -to-noisa  ratio  is 


For  regular  cross-correlation  detection,  the  second  term  is  not  there  and 
tha  Rg(L)  is  not  there,  leaving  2NoRs(0)/T.  For  regular  auto-correlation, 
corresponding  to  L  *  0,  wa  have 

4WN?/T  f  4N  (0) /T  ns  the  value  for  u*. 

O  O  ■  K 
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: - sin  C4"  TOT 

*  h  'i*wr; — 


t  Rs(L)/Rs(0)* 


Actually  when  G(f)  <<  Nq  -  say  by  at  least  several  db?  e.g.,  when  (S/N^ 

is  perhaps  -lOdb  -  the  last  noise  term  isn’t  important  either.  Then 


and 


ainUf'WL)  ’ 
MT — 


ls\  .  _ (s\2  r,  .  ain(4irWL)V1 

Wo  ■  2™Wi  l1  -Twr-J 

Maximum  (S/N)q  occurs  for  L  =  ( 3/8W) ( 257 . 5/270) .  This  is  4  db  better  than  the 
value  for  L  ~  0.  With  L  -  (3/8W) ,  corresponding  to  an  angle  of  270  degrees 
in  the  sine  function  above,  there  is  negligible  reduction  in  (S/N)o  from  the 
maximum.  We  note  that  for  L  equal  to  an  integer  multiple  of  (1/4W)  and  for 
L  -  oo,  the  improvement  in  signal -to-noise  ratio  over  the  value  for  L  -  0  is 
3  db. 
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